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Abstract 

We estimate the size of the triple Pomeron vertex in perturbative QCD and compare 
with the phenomenological value extracted from Regge fits to experimental data. For 
simplicity, the results of the QCD analysis are taken in the large- A'^c limit. We find 
that the perturbative triple Pomeron coupling is of the same order of magnitude as the 
observed one. We also estimate the size of the Pomeron self energy and its contribution 
to the renormalization of the Pomeron intercept. The effect is very small, in agreement 
with previous nonperturbative estimates. 



1 Introduction 

The Regge description of hadronic high energy scattering processes contains a few funda- 
mental parameters which are of nonperturbative nature. Their values have been extracted 
from the analysis of a large variety of experimental data, and, so far, there exist no calcu- 
lations within QCD which would allow a comparison of theory and experiment. Prominent 
examples are the Pomeron intercept ap{0) ~ 1.08 and the Pomeron slope a'p ~ 0.25 GeV~^, 
seen in the total cross section and in elastic scattering, and the triple Pomeron coupling g^p, 
defined and measured in high mass diffraction. 
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Whereas the former two parameters refer to the (effective) Pomeron seen at present en- 
ergies, it is widely beheved that the latter one provides information on the origin of the 
Pomeron: a Pomeron with intercept exactly at one would exhibit features that are typical 
for systems near a phase transition point @, ^ . In such a situation the triple Pomeron vertex 
which describes the splitting of a single Pomeron into two Pomerons then provides the start- 
ing point for calculating correlation functions, critical indices etc. For example, using a field 
theoretic description of the Pomeron, the value of the triple Pomeron coupling determines 
the size of the self energy, it renomalization of the intercept etc. In reality, the intercept is 
close to unity (but not exactly at one), so it is likely that, at present day energies, we are in 
the vicinity of a phase transition, and the triple Pomeron vertex plays a fundamental role. 



In perturbative QCD, the Pomeron is approximated by the BFKL calculation |1|] (in LO and, 
more recently, also in NLO 0). However, the values for Pomeron intercept and slope are 
not very close to the observed hadronic values; moreover, the BFKL approximation can be 
justified only for scattering processes in which the scattering objects have a small transverse 
extension (7* — 7* scattering, or onium-onium scattering). As to the Pomeron slope, for 
t ^ 0, the BFKL amplitude predicts a small value, whereas at t = the t-slope is singular, 
reflecting thus the long distance behavior of the perturbative massless gluons. The next pa- 
rameter, the perturbative triple Pomeron vertex, has first been calculated in P, ^, starting 
form the high energy behavior of QCD Feynman diagrams. Later on, independent deriva- 
tions have been performed, within Feynman diagrams |^ ^, using a Wilson line approach 
and within the QCD dipole approach [0, |ll[]. As far as the numerical computation 



of this perturbative coupling and its comparison with the experimental hadronic vertex is 



concerned, an important step has been done in |jT2|, |T^: the analytic expression derived from 
the underlying Feynman diagrams contains conformal integrals which have been computed in 
||T^ , p!3| . These results, however, do not yet allow for a direct comparison with experimental 
data: as it was the case already for the BFKL approximation in elastic scattering at t = 0, 
also the perturbative triple Pomeron vertex has a singularity at zero momentum transfer. 
Any numerical estimate, therefore, will depend upon the way in which this singular behavior 
is treated. 



There is no doubt that perturbative QCD cannot be used in hadron hadron small-angle 
scattering. Nevertheless, the analysis of perturbation theory in this high energy limit pro- 
vides the first step towards the 'real' theory, and it is important to see, 'how far away from 
reality' we are in pQCD. It is the purpose of this paper, to attempt a numerical estimate of 
the perturbative triple Pomeron vertex and to compare with the hadronic value. We start 
from the Feynman diagram analysis of 0, 0, and we make use of the numerical values 
obtained in |jT2|, |13|. For reference we use the cross section formula for diffraction in the 
triple Regge region: we derive a value for the perturbative triple Pomeron vertex which 
can be compared to the measured hadronic value. We also estimate the self energy of the 
BFKL Pomeron. Some of our results differ from earlier estimates, contained in the literature 

ra, PI, HI. 
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When trying to compare BFKL predictions with the Pomeron parameters measured in 
hadron hadron scattering, we will face a few difficulties of general nature. First, hadron 
hadron scattering, to a very good approximation, has been parametrized by a simple Regge 
pole in the complex angular momentum plane; the leading BFKL singularity, on the other 
hand, is a fixed cut which leads, in addition to the Regge exponents to logarithms of the 
energy. Furthermore, BFKL scattering amplitudes are slightly singular when the momentum 
transfer t is taken to zero; this singularity reflects the 1 / ^^-singularity of the zero mass gluon 
propagator in perturbative QCD. Nonpcrturbativc effects, therefore, are expected to be par- 
ticularly strong near t = 0, and a comparison between perturbative Pomeron parameters and 
the measured nonperturbative values looks more promising in the region of nonzero values. 

Our paper is organized as follows. In section 2 we briefly review, for comparison, the pertur- 
bative BFKL Pomeron in elastic scattering. In section 3 we turn to the triple Regge region 
of diffraction and define what we mean by a 'triple Pomeron vertex' in perturbative QCD. 
Section 4 deals with the self-energy of the BFKL Pomeron. The numerical evaluation will 
be done in section 5. In the final section we give a summary and a few general comments. 
Some technical details are put into two small appendices. 



2 Elastic scattering 



In order to find the correct normalization of the triple Pomeron vertex we have to start from 
elastic 2 — > 2 scattering. Let us write down the Regge ansatz for an elastic 2 — > 2 scattering 
process. For a Pomeron pole in the complex angular momentum plane the elastic amplitude 
is 



i-ei = -e - - >g^s ( — ) , (1) 



the elastic cross section has the form 



da 1 4 /sX2a(t)-2 

9n{-) , (2) 



dt IGtt vsq 
and the total cross section is 

/ „ X a(0)-l 

In the following we want to compare these expressions with the ones obtained from 
perturbative QCD in the Regge limit. In particular, we have to relate the residue functions 
gijv to impact factors which naturally arise in a perturbative analysis. For reasons which will 
become clear soon, we will have to define forward and non-forward coupling functions, gp 
and Gnf, resp. 

For a 2 — > 2 scattering process (e.g. gluon-gluon scattering or 7* — 7* scattering) a color 
singlet exchange leads, in lowest order ag, to the form 
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All the momenta are living in the transverse plane, q (with t = —q^) denote the momentum 
transfer and s the squared center of mass energy, resp, and is the impact factor of the 
scattering particle i. As an example, with this choice for the integration measure, the gluon 



impact factor can be written as g'^2y/Ti Nd ^ N"^ — 1. 

Summing all the contribution in the leading logs approximation leads to the BFKL 
Pomeron exchange; instead of the two gluon propagators we insert the BFKL Greens func- 
tion. The new amplitude reads 

A^^t = f / (03^'^i(^' 9 - ^)G{y\k, q-k;k\q- k')Uk\ q - k'), (5) 

where y is the rapidity variable. Clearly, for — ^ 0, when all rungs of the BFKL resumma- 
tion decouple, this expression reduces to the two gluon exchange, which means that we use 
the following normalization: 

\im G{y\k, q-k;k',q- k') = -J^^^5^^){k - k') (6) 
ois^o k[q — kY 




(a) (b) 



Figure 1: Elastic process 



2.1 The BFKL Pomeron Green function 

Eigenfunctions of the BFKL kernel, E^ j^, are well known in coordinate space, where its form 
is dictated by conformal invariance 



h 



10' 20, 



with rio = ri- tq etc, h = {l + n)/2 + iv, h = (1 - n)/2 + iu {h* = 1 - h, h* = 1 - h), 
and standard complex notation for the two-dimensional vector is used on the right-hand 
side. Fourier transforming (we use the Lipatov's convention which assigns a l/(27r)^ to any 
coordinate integration) to momentum space one finds |TH 



E,n{k^,k,) = J Sj^^/.,;.(^i>^2)e^(^^-^^+^-^^) = C{x{k„k,) + {-irX{k„k,] 



(27r)2 (27r) 
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The coefficient C is given by 

C = ^j-^hh{l-h){l-h)T{l-h)T{l-h). (9) 
The functions X in complex notation can be expressed in terms of hypergeometric functions: 

X(...M^(|)'"(|)"V(l-.,2-/,2;-|)F(l-/,.2-M;-g • (10) 

This analytic form does not contain any term of the type 5^(fci) or 6'^{k2) which are present 
in the coordinate representation (|^). For the impact factor of a colorless external particle we 
have the well known property, that it vanishes for zero gluon momentum: in this case the 
delta-function type contributions do not contribute. For simplicity, we therefore will ignore 
them. 

The Pomeron intercept has the form ap{0) = 1 + x(z/, n) where 

x{v,n) = aAmi)-ij{^^+iu)~ij{^y-iv)\, a, = (11) 



TT 



Let us now consider the Pomeron Green function in the coordinate representation 
Making use of the Casimir operator properties of the Mobius group, one can chose a rep- 
resentation which is more convenient to perform a Fourier Transform. In particular, after 
considering the relation 

= 16 y + (n - l)V4][z.2 + (n + 1)74] WWE.^j^ip.o. P2,){12) 

one can write 

Gi^\y\pi,P2;pvP2') = J rf/irfVoe^^('^'"^Ar,|aina2pi?,,S,(pio,P2o)i5^u(pi'o,P2'o) (13) 

where we use the measure in the conformal weight space J dp = J2n I dv with the following 
normalization factor Nh. 

AT, = + (14) 

[z/2 + (n- 1)2/4] [z/2 + (n + 1)2/4]' ^ ' 

In order to find the momentum representation we take the Fourier transform taking into 
account the total momentum conservation, and we obtain: 

Gt\y\k,M;kvk2:) = (271^1 rf/ie^>^('^'")A^;,x(27r)7fcinq-fci|2E,,^(fci,g-fcO%,(fci',q-fci' 

(15) 

The (27r)^ factor in front of the integral comes from the normalization (P); the second (27r)2 
factor results from the Pq integration, together with a 6^'^\ki -\- k2 — ky — ^2') related to 
the overall momentum conservation. As before, q = ki + k2 = ky + is the conserved 
exchanged momentum. This form of the Green function is amputated on the Ihs, i.e. for 
the gluons with momenta fci and k2. Clearly, dividing by |fci|2|qr — fcip one arrives at the 
nonamputated Green's function 6*2^"^''. 
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2.2 Extraction of the couplings gp and gjyp 

Substituting in the expression of the BFKL Pomeron Green's function the elastic scat- 
tering amphtude can be written as: 

A^i'iq) = f (2vr)^E / rf^Ar,e^^(^'")$^(g)$^*(g) , (16) 

where 

^^liq) = J ^jsEkAk, q - k)Mk, q - k) (17) 

are the impact factors in the conformal representation, i.e. integrated with the BFKL 
Pomeron eigenstates, and q is the total transverse momenta exchanged. 

For our purposes it will be sufficient to consider the elastic scattering of identical particles. 
We will consider the forward {q = 0) and the non forward {q ^ 0) case separately. Since we 
are interested in the leading high energy behavior, we restrict ourselves to the conformal spin 
n = 0, and we perform the integration in u in the saddle point approximation for y oo. 
We therefore need to know the behaviour of "I'f as a function of u. In the appendix we show 
that the forward ($f) and non forward {^nf) cases are different. Near the saddle point at 
u = we find: 

$F = —%F + 0(z/°) , <l>^^ = %NF + 0{u) , (18) 

Moreover, we have = IQu^. Next we need the expansion 

X(i',0) = Xo-ai^\ Xo = 41n(2)a, a = 14C(3)a,. (19) 
In the forward case g = we obtain 



^i^ = -8(2-)^|^^oi.e-o, (20) 
while in the non-forward case we have 

A^^ = ^^8(2^)'|^^'^oiVFe^"" • (21) 

The transition between the non-forward and the forward region is a delicate matter. The 
different large-y behavior {y~^^'^ and y~'^^'^) of the Ap and A^p amplitudes originates from 
the different small-i/ behavior (note, in particular, the 1/v singularity in the forward impact 
factor (18)). As we will show in the Appendix B, this 1/v singularity comes from the large-ro 
domain, and it reflects the perturbative nature of the BFKL Pomeron. In particular, it is 
related to the singularity of the perturbative gluon propagator at zero momentum, and it 
must disappear after the introduction of an appropriate infrared cutoff Considering the 
elastic scattering amplitude as a function of the momentum transfer t, the difference in the 
small- z/ behavior of the non- forward and the forward results leads to a cusp at t = 0: the 

■^This was demonstrate using a simplified form of the BFKL kernel in 
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elastic cross section (2) has a finite hmit ed, t — 0, but its ^-derivative at t = is infinite. 
Generally speaking, in the BFKL approximation the point t = exhibits the perturbative 
nature most explicitly, and changes from perturbative to nonperturbative QCD are expected 
to be most dramatic in this kinematic region. 

After these general remarks we are now able to extract, by comparing (20) and (21) with 
(1), the Regge residue factors gp and gNF' 

a^p = 2'"%„p(2^r" (j^Y (22) 

The fact that these couplings have a residual y-dependence is a consequence of the branch 
cut nature of the BFKL singularity in the angular momentum plane: eqs.(l) -(3) are valid 
for Regge poles. We shall use these relations in the next section in order to extract the triple 
Pomeron vertex g^p. 



3 Triple Pomeron amplitude 

Having collected all necessary ingredients we now turn to the central topic of this study, the 
triple Pomeron vertex. We again start from the Regge form for the diffractive cross section 
in the triple Regge region, assuming Regge pole singularities in all three exchange channels. 
The cross section is obtained from the 6-point amplitude (Fig. 2) by taking the discontinuity 
in the diffractive mass squared, M^. We define 

Regge theory gives 



1 



-ifa(t)|2 n"^ n ( 




9f9%f9^p[th] — (24) 



where g^p is the triple Pomeron vertex, t denotes the momentum transfer, and sq is an 
energy scale. The triple Pomeron vertex depends upon t. It will also be convenient to intro- 
duce the rapidity variable Y = log(s/so) and the rapidity interval of the diffractive states 
Ym = log(M^/so)- It will be our aim to extract the counterpart of g^p in the framework 
of perturbative QCD (by analyzing, in the leading log approximation, the analogous high 
energy limit in perturbative QCD), and to compare its value with the empirical value in pp 
scattering. To this end we consider a (hypothetical) process in the triple Regge region, e.g. 
7*7* — * X'-f*, for which the use of perturbative QCD can be justified. Because of Regge 
factorization, the value of the triple Pomeron vertex will be independent of the external 
particles (e.g. 7* with virtuality Q^). 
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3.1 The QCD amplitude 



Let us look at the analysis of QCD Feynman diagrams in the leading log s approximation 
and recapitulate the main results. To this end we write the general integral representation 
for the differential cross section of the 3^3 process in Fig. 2 (see eq. (2.2) in 0) 

^.oda^'^'ff^ 1 r doj f dui f du2 ( M^Y f s \^i+^2 ^ x . ^ 

*'W = Ife/2^l/d/dM U) «".e,fK-.,-..0,M), (25) 

where the are signature factors. 




Figure 2: Diffractive amplitude 



The differential cross section (^5]) can be represented by diagrams built from reggeized 
gluons. The set of diagrams which is of particular interest for us is illustrated in Fig. 3a: el- 
lipses denote impact factors, circles the two-gluon BFKL Green's function, and the triangle 
the 2 — > 4 gluon vertex The shaded box denotes the evolution of the 4 gluon state, medi- 
ated by the sum over all pairwise interactions between the four reggeized gluons. The last 
interaction (in Fig. 3a: the lowest one) has to connect one of the two reggeized gluons on the 
Ihs with one of the gluons on the rhs. Let us briefly recapitulate how this result is obtained. 
The analysis of Feynman diagrams in the high energy limit leads to gluon amplitudes D2-, 
-D3, and which satisfy a set of coupled integral equations (Fig.4). These functions are 
nonamputated, i.e. they contain reggeon denominators for the outgoing (reggeized) gluon 
states. Removal (amputation) of these reggeon denominators leads to the corresponding 
functions Ci [i = 2, 3,4). For example, 



D 



(26) 



In order to obtain the partial wave F of the triple Regge cross section, we attach two 2^2 
BFKL Green's functions to the amputated function C4, one for the two outgoing gluons on 
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the Ihs, another one for the two gluons on the rhs. In order to avoid double counting we have 
to require that the last interaction inside the four-gluon state has to connect one of the two 
gluons on the Ihs with one of the gluons on the rhs. As a result, we arrive at the following 
expression for the partial wave F in the triple Regge cross section formula: 



cu-f3i'- P2' u;-(3^-(32-/33'- Pa 



)G2^2il2,CO,)G2^2iM,U2) 



(27) 

where (3i are the trajectories of the reggeized gluons, ¥2^2 is the BFKL kernel (without the 
gluon trajectory function, but including its tensor color structure, here acting for the pairs 
(12) and (34) which are in a color singlet state), and (^2^2 denotes the full non forward 
BFKL Green function. A convenient way to rewrite F in terms of D4 is: 






Figure 3: Triple interaction 



D 



[UJ 



H^) + Df^ ® (r(13) + V{IA) + r(23) + V(2A)) 



5^2^2(12, a;i)G2^2 (34, u;2) 



(28) 



where H/^ is the standard BKP evolution operator for the 4-gluon state. 

As the last step, one has to apply a reduction procedure to the amplitude 1)4. As a part 
of the coupled integral equations, still contains reggeizing pieces: the outgoing four gluon 
state may contain configurations where a pair of two gluons is in an antisymmetric color 
octet configuration, which satisfies the BFKL bootstrap condition and collapses into a single 
gluon. It is convenient to remove these configurations, i.e. to define amplitudes D^^ which 
are irreducible with respect to the bootstrap property. This reduction has been described in 
0, and decomposes into the two terms = D^ + D{, separating the reggeizing (R) and 
irreducible (I) parts. In (p8|), let us first consider the irreducible part, D[. As shown in [Q, 
D\ consists of the diagrams of Fig. 3a which we have described before. The triangle - with 
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two gluons entering from above and four gluons leaving below - defines the triple Pomeron 
vertex, and its structure is quite simple: 

Sbb' {5a,aJa,a,V {12, 34) + Sa.a^Sa.a.V {13, 24) + 5a,aA2a,V{U, 23)) , (29) 

where the b, b' are the color labels of the reggeized gluons of the ladder above the triple 
Pomeron vertex, the color indices of the reggeized gluons inside the two lower ladders 
(counting from left to right), and the arguments of the function V refer to the momenta 
of the gluons. Below this vertex, before the two gluons on the Ihs and the two gluons on 
the rhs are restricted to color singlet states and branch into the two disjoint BFKL Green's 
functions, all pairwise interactions between the four gluons have to be summed. However, it 
is easy to see that any rung between two color singlet two-gluon states costs a suppression 
factor of the order 1/N^'- in the large- iVc limit, therefore, in Fig. 3a the interaction inside the 
shaded area can be neglected, and we are left with the diagrams of Fig. 3b. In (|28|) , the factor 
uj — Hi cancels the evolution inside Df, and the terms proportional to ^(13) etc. drop out. 
In (|29|) , only the first term contributes to the large- iVc limit. 




Figure 4: Chained equations for the multi reggeized gluon amplitudes. 

Before we write down the explicit expression for these diagrams, a few words about the 
contribution of Df. It is convenient to go back to (^) and (p7|). Since Df is nothing but 
a BFKL ladder in which, at the lower end, the reggeized gluons split into two (or three) 
elementary gluons, it provides an extra contribution to the triple Pomeron vertex. From 
the color structure of D^, given in 0, eq.(4.3), it can be shown that this contribution is 
subleading in 1/Nc. In conclusion, the large- iVc limit therefore reduces the diffractive cross 
section, given in 0, to the diagrams shown in Fig. 3b which is very similar to the 'Regge 
pole' diagrams of Fig.3c. 

3.2 Extraction of the triple Pomeron vertex 

Let us write down the analytic expression for the the generalized six point amplitude a^. 
oe = '^sj (27r^3^ ^{2n)^^' ^^^^^^^ ~ ~ ^2')'^c.(q, L - q)Gi'^\YM\q, L - q, Qi', ^2') 

J (2n)^ -ki- k2) o"- '{l^ - fcs - k4}V2^4{qi,,q2'\ki,k2;k3,ki) 
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I ^^2^^^^^^ ~ ^"''1^'' ^2'' ~ ^3')'^'7(fc3', ^7 - ^3'), (30) 

where Zj (with — tj) are the momentum transfers through the impacts factor $j. Later, on 
we will take to, = and define the triple Regge cross section. 

Following the procedure of the previous section, we perform the Fourier transform for the 
ta-channel (last two lines in (|30|)); we write a representation of 5^'^\la — Qy — ^2') together 
with the Green function in the last line, see(|T3p, to obtain 

- Qv - Q2') &2'\yM\q, la - q, Qv, 92') = {^^?J rf/ioe"'^^^("°) iV,„ x 

/rfXe-^"-(^"-9i'-^^')|gi,nq2'r4.,;..(q,/a-q)^L,;i.(9i',g2')-(31) 

The other two t — channels with their (non amputated) Green functions and 6 distribution 
can be rewritten in a similar way (using an expansion in the conformal weights /i/j and hj); 
the only difference is the absence of the factors |A;i/p|Z^ — fci/p and |fc3/p|Z^ — fca/p. Using 
the conformal representation of the impact factors $f'(Zj) for i = 7 given in (0), we 
obtain the following expression for 



2s J dfi^dfipdfi,e''^'^^''-^+^^-''^'^^^^^^^^^ n d'P. 

1=0, /3, 7 

e^^P^-^^^P^-^^-P^-^-^ J nci^%4„;i,(fei,fe2)e-^^/^-(^^-^^^)4„;i,(fc3,fc4)e-^^v(fc3+fc 
J d'qi,d'q,, V2^Mv, q2'\ki, k^, k,, k,) \qy\^\q^,\^ El^^j^^{qy , q,,)e'P^'<1^'+1^'^'^ (32) 



Before we evaluate the Fourier transform of the last line in (|5^, we note several simpli- 
fications. First, the 2 — > 4 vertex will be simplified by the fact that the gluons (1, 2) and 
(3, 4) couple to two BFKL pomerons in color singlet states. This fact considerably reduces 
the number of contributions coming from ¥2^4. only four identical contributions are left. 
Next, we restrict ourselves to the large limit which eliminates the nonplanar part of the 
2 — »• 4 vertex. As a result of these simplifications we can write: 

J d^qyd^q2,V2^Mv^q2'\ki,k2\k'i,ki) \qy\^\q^,\^ El^j^{qy , q^,) = 

Civ J d^qyd^q^, A2-,^{ky fcg + fcs, ^4191', ^2') l^i' 1^192' P^/I,/,,(9im ^2') (33) 

where the constant 

r< A_ ^4 ^ ^4 ^^c ^ o3_7/2^2 tQA\ 

32 /N2 _ I 2 /N2 _ 1 
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follows from our choice of the normalizations we have made for the integration measure and 
for the impact factors. Details are given in appendix A. As to the color factors, we keep only 
the leading term in the large Nc limit. 

After performing the Fourier transform of the last line in ( p2D 0, |^ (i.e. performing the 
transition ki pj) one obtains for the triple Pomeron coupling 



1 


PlA 


2 


|Pl2 


2 


P24 


2 



(2vr) 

^5^^){p^^)lQKK{l - K){1 - K) ^ A ,hM-^P^o.) (35) 

(27r)* |pi2r|P24r|P4ir 

The Fourier transform of the remaining fcj dependent part is easily done and leads to the 
two factors E^^-^^ipip, P2/3) and E^^-^^ips^^p^^). 
As a result, we can write 



7 X 



/ I ^21 ^121 ^]2 ^hf,,hAPlf3, P2f3)Eh h {p2-r, Pi-f) E^^^h^ {pia, Pla), (36) 



where C2V = Ciy/(2vr)'^. 

The integral in the last line of (^) has been calculated explicitly in [|T^, |r^, where the 
conformal invariance has been used explicitly. The result can be written in the form 

Q{1 - ha, hp, h,) {p-J-'^'p-^-^p-p^'^ X {h.c.)) (37) 



where the function n can be found in |T^, |12|. The exponents are defined for the general 
conformal covariant three point function: Aij = hi + hj — hk^ij and remembering to use for 
the index the weight (1 — ha) = h*^, which is due to the fact that one function is complex 
conjugated. That means A^^ = 1 — ha + hp — h^, etc. 

We shall now consider the limit Zq, = (keeping tp = — and = — still independent 
from each other) and perform the remaining p integrals still present in (pB]). Explicitly we 
have to calculate: 

/, = / d'pj'ppd'p^ (p^^V^f-p^,''^^ X {h.c.)) e^PM^P^l-' (38) 
The integration over p^ can be done easily: 

' d'pa (p:^V^,^"^ X {h.c.)) = f{ha, hp, h,) (p;-^"^-^- X {h.c.)) (39) 

^Comparing with an analogous expression in jl^ , we differ in the Casimir operators which are not present 
in O 
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where 

f(h h h ] - ( 1^A^.-A„, ^(l-A^^)^(l-A^^)^(l-Aaff)^(l-A^^)sin(7^A»;3)sin(7^Aa^) 
^' ^ ^ r(2-A,^-A„,) r(2-A„^-A„^) sm(7r(A,^ + A,,)) 

- ( 1 )^(h-y-hf,) -hf} + h^)T{K -h^ + hp) T{K -hp + h^)T{ha - + hp) 

^ ' T{2h^) V{2h^) 

sin(7r(/ta -hf) + h^)) sin(7r(/ia - h^ + hp)) ^^^^ 
s\.n{2Tihct)) 

Therefore one is left with the following integration 

Performing a shift in one of the integration variables, one integration can be done to extract 
the momentum conserving 5 function: 

h = {2'K)H^^\lp + L){-l)^^-<+~^^'<f{h^,hp,h^) j rf2^ X {h.c.))e'P^^. (42) 

Also the remaining integral can be done easily. Introducing ^ = 1 — Aq,^ — A^^ ~ ^p-y and 
^ = 1 - A„/3 - A„^ - A/3^, we have 

/ d'pp^p*^e'P^^ = 2n 2i+«+f ^^WZ^ (43) 
Returning to (pHDwe have 



h = i2n)'6^'\h + l.)fih^, hp, /.,)(-l)^^^+^^n«-« 2i+«+f (44) 

with 

^ = I - {{I - ha) + hp - h^) - ((1 - ha) + h^ - hp) - {hp + h^ - {1 - ha)) = ha - hp - h^ 

(45) 

and 

^ =ha — hp — h^ ^* = —I — ha + hp + hj C = — ha + hp + h^ (46) 
Collecting all this in the formula for the amplitude in (p6D , we obtain the result: 



a6= 2sC^v j dpadppdp^e^-'^^''-^+^'''''^'^^-^^''^^+^^^^^^ x 
16haha{^ — ha){l — ha)^{l — ha, hp, hy) f{ha, hp, hy) X 

A,,+A,,^e--c 2i+c+e ^^l+J) (47) 

Here we have removed the 6^'^\lp + ly) function. The constant C^v has the form 

Csv = {27r)'C,v = ^ (48) 
13 



Let us now consider the saddle point approximation, assuming that both Ym and Y — Ym 
are large. Following the standard BFKL arguments, the leading contributions will come 
from the conformal weights being close the value 1/2, i.e. from conformal spin equal to zero 
and small i^'s. In this region we have 

hi = ]^ + iui, ^ = -1 + - - z^-y), (49) 

and therefore 

. _i X X , . , (50, 

1(1) 1(1) —smlmua 2,Va 



(--LJ ^ */3 2f(|)^ ^ 

Moreover, 16/iQ(l — ha)ha{^ — ha) ~ 1 and ~ 16z/^. We note that the function /, near 
la = 0, behaves in z/^ in the same way as the forward impact factors (18) Therefore, with 
respect to the ta channel we are facing the same problem as discussed in the previous section. 
Substituting this behaviour in (|i7| ) one has 



2'2'2y 4I/J ^' 



where 



^hfl*^h-.,*^l^^-l-2i{uc-ui3-uj) ^52) 

CUv ^ 2",n (1. i, 1) C,y ^ 2Hl (i. 1,1) ^ (i i i) . (53) 



The numerical value for fi(|, |, |) has been found in [|T3], |T^ : 

(54) 

The remaining part of our saddle point analysis of is analogous to what has been 
done in subsection 2.1: the integration over i>a goes with an impact factor in the forward 
direction, whereas for the two integrations in i/^ and we have the freedom to vary tp 
and tj. following our discussion after (21) we chose to stay away from the 'most dangerous' 
points tjs = = 0, i.e. we perform our comparison in the 'safer' region tp = ^ 0. Putting 
1/3 = la = I and using (18) we find that the dependence in I in the dominant contribution. 
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as selected by the saddle points at t'a = 1^/3 = 2^7 = 0, is just The integrals are trivially 

done and lead to the result: 



Using the relations in (^) for rewriting the impact factors in terms of gp and gNF-, we arrive 
at: 



1 



a6~^^/F^7;vF 27/2 (2vr)i5/2 UanJ 1^ [2a(y - n,)] V ' ' ^^^^ 

Finally we have to relate the generalized amplitude ag to the triple Regge cross section 
formula: 

(ia^''-^ff^ 1 

Comparing our result with (p^, we are able to extract the 



g^p 



C4V 1 / 27r 




|/| 25/2 (27r) 15/2 V2aFM. 

L 2^ /I 1 1\ / 27r \* / 2n 



|/|(2vr)^"^^^^l2'2'2jl2fl^j I, [2a(r - Fm)P j " ^''^ 

We note that this expression is proportional to cel^'^, i.e. there is a very mild dependence on 
the strong coupling. Considering ~ 0.3 we get 

The strong dependence upon the momentum transfer near 1 = which is closely connected 
with the perturbative zero mass gluon confirms our expectation that a comparison between 
perturbative and nonperturbative Regge parameters can be done only in the region of finite 
momentum transfer, and one cannot expect more than an order-of-magnitude estimate. 

Let us comment on other results of this vertex contained in the literature. Within the 
dipole picture the triple Pomeron has been derived in [0: an expression for the triple 
Pomeron coupling can be derived from eq.(61), but no explicit expression or numerical 
number Vq has been given in this paper. In explicit expressions for the triple Pomeron 
vertex can be found: our result disagrees, both in the energy dependence and in the overall 
normalization. The result of [|1^ is closest to ours, but, again, we disagree in the overall 
normalization and in the energy dependence. 



4 The Pomeron self energy 

As an important application of the triple Pomeron coupling we estimate the size of the 
Pomeron self energy inside a 2 — > 2 scattering amplitude. To this end we replace, in the 
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amplitude of the previous section, the two impact factors and by another triple 
Pomeron vertex V2^4 which though a BFKL Pomeron and another impact factor couples 
to the lower external particle. As a modification of Og, one has to consider an additional 
momentum integration over the Pomeron loop, in the momentum variable Ip = la — l-y- For 
simplicity we will restrict ourselves to the forward direction = 0. Denoting by ^' the same 
combination (39) of conformal weights as ^, with ha being replaced by ha', the lf3 integral 
can be carried out and leads to the conservation of the conformal weights above and below 
the Pomeron loop: 



j2j I j--ha+hp+h.,~lj*-ha+hi3+h^-l\ ( j* h^, -hp-h-, jh^, -hp-h-, 
d 1/3 Uf^ Ip I Up If^ 

j2i jh ,-ha-lj*h i-ha-1 (271") r r/ \ / p n\ 

d Iplif " = -^On^n^,d[Va - Ua') ■ (60) 

Another important ingredient to the Pomeron self energy is the minus sign relative to the 
BFKL amplitude (5). 

We want to evaluate this loop correction to the elastic scattering amplitude, using again 
the saddle point approximation. Starting from the expression (|52D, inserting the result ( |60D 
and performing the integration, we can write 

2 4(27r)4 \16{27ifjJ 

noting that the z/^ and z/^ integrations give identical factors. Let us also note that all the 
rapidity intervals must be large enough to allow the application of the leading log approx- 



imation. This expression can be compared with a similar result in The remaining u 

integrations give: 



XO {Y1+Y2) 




2x0 (y-Yi-la) (an) 

[2a{Y - n - Y,)f 

This result represents the one loop self energy correction to the BFKL approximation 
(pop. Apart from the fractional powers of ag in front of the rapidity factors, the overall power 
of as inside the C4V'-factors is a^: compared to the LL BFKL approximation our expression 
is down by two powers of a^, i.e. the self energy correction belongs to NNLO and thus is 
beyond the NLO corrections calculated recently. Note, however, the exponent 2xo in the 
last line: for large rapidity intervals Y — Yi — Y2 this energy factor renders the one-loop self 
energy correction more important than the NLO corrections to the BFKL kernel. 
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5 Numerical estimates 



5.1 Triple Pomeron vertex 

In this final part of our study we use our analytical formulae to obtain numerical estimates. 
We begin with the phenomenological vertex extracted from the pp p + X data in 
the framework of the old triple-Regge analysis |21|, ^ and compare with the perturbative 
triple-Pomeron vertex g^p. 

It was observed that the t-dependence of the triple-Pomeron contribution to the diffrac- 
tive dissociation cross section is consistent with the t-behaviour of the proton-Pomeron vertex 
square; that is M'^da/dtdM'^ oc g%{t). Hence the t= -F-dependence of the triple-Pomeron 
vertex gspit) must be small; for the t-slope of g^p we estimate B^p < 1 GeV^. Therefore 
we may consider relatively large |I| ~ 1 — 2 GeV where, in the perturbative calculation, we 
are away from the QCD dangerous region 1 = 0. In our normalization the phenomenological 
analysis gives ^ 

gsP ^ 0.5 - 1 GeV'K (63) 

Note that when applying our formula (23) to experimental data and extracting a numerical 
value for the triple-Pomeron vertex, this vertex has to be viewed as an effective vertex, i.e. 
it already accounts for screening corrections due to multi-Pomeron cuts. So the bare vertex 
may be larger by a factor of about up to 2 - 4 Thus, at the experiment we "observe" a 
bare vertex of the order g^p ~ 2GeV^^. This value corresponds to the events with a gap 
size AY = ln{s/M'^) = Y -Ym between 3 and 5. 

Turning to the perturbative analysis, we first note that, in order to justify the saddle point 
evaluation of our integrals, we need AY > 4. This is just the region of z = yasNc/71 > 1 
where the asymptotic component of the BFKL solution (with conformal spin n = 0) starts 
to exceed the lowest order two-gluon exchange contribution. At the same time the width of 
the saddle point Su ~ l/(aAy)^ ~ 0.3 (for = 0.3) becomes sufficiently small. Thus, if we 
choose Ym = AY = Y — Ym = 4, = 0.3 and 1=1 GeV~^, we obtain g^p ~ 0.6Ge\^~^. 
It follows from (59) and from our discussion before that this value has large theoretical 
uncertainties: changes in 1, Y and Ym have a stronger influence on the numerical value of 
the perturbative triple Pomeron vertex than on the nonperturbative triple Pomeron coupling. 
As the result of our analysis, we therefore present the range 

g3P ~ 0.2 - 1.7 GeV'^ (64) 

which is related to the ranges of values 3 < AF < 5 and 0.5 Gel^"^ < 1 < 2GeV-\ 
Surprisingly, these numerical value are not far from the experimental value discussed before. 

Finally, we would like to mention that this perturbative value may be a little overes- 
timated. Namely, we have to remember that, in order to arrive at the triple Regge cross 
section formula, we had to take the discontinuity of at = 0. In analogy with the impact 
factor (18), the triple Pomeron vertex at has a singularity ~ 1/z/ just at the saddle point 
z/ = 0; as we have discussed in section 2, this singularity is a result of the massless gluon 



^In particular, the gap survival probability within the ISR energy domain, calculated in the formalism of 
ref. pi, is equal to = 0.25 - 0.33. 
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propagator, and it should disappear after the introduction of an appropriate infrared cutoff. 
Numerically, the presence of the 1/v singularity should lead to enhancement, and our value 
of the perturbative triple Pomeron vertex may, in fact, therefore be overestimated. 

It is interesting to compare our result with the numerical studies of the Balitsky-Kovchegov 
equation ^ (BK) equation. The equation includes the LO BFKL evolution and 

accounts for the triple Pomeron coupling summing up the fan diagrams in terms of the 



dipole-dipole interaction. In the recent paper the begining of saturation was observed at 



rather small F ~ 2 - 3. However first the saturation is reached for a large size dipole where 
the absorptive corrections are much stronger. This is a dangerous region. Even without the 
confinement and for a fixed as coupling we faced here two problems: 



on one hand, at small Z ( A; in the notations of ||26|) the l/u singularity plays a crucial role , 
as it was discussed above, 

on other hand, the whole approach can be justified only for the case when the rapidity inter- 
val Ay occupied by each Pomeron is large enough, there is no this condition in BK-equation 
and the Pomerons can split immediately, especially for a large size dipoles. 

To avoid these problems we focus on dipoles of a smaller size, smaller than the initial 



(input) size l//co, taken in to be 1 GeV~ . Here the absorptive effects reveal itself at 
y ~ 4 - 6 (see Fig.2,4 of ref.||26|) . This is in agreement with our expectation. Based on 
the simplified form of the first (order of (?3p) fan diagram contribution (23,24) and taken 
the efective perturbative vertex Q'^p ~ 0.6 GeV^^ (which was evaluated just for AY ~ 4), 
we found that the triple Pomeron amplitude becomes comparable with the single Pomeron 
exchange at y ~ 4 - 5; we choose Z ~ 1 GeV, Q]^ = 10 GeV~^ (corresponding to 0"*°* = 40 
mb and = 0.2 (corresponding to the LO BFKL ujq = 0.56) as it was done in E^). 



5.2 Renormalization of the Pomeron intercept due to the Pomeron 
selfenergy 

Finally, we estimate the size of the Pomeron self energy correction and its influence on the 
intercept of the perturbative BFKL amplitude. In the loop amplitude (|62D we still have the 
integrals over the rapidities Yi, Y2, and the dominant contribution comes from the region of 
small Yi, Y2, where AY = |y — Yi — 1^2! Y. This limit corresponds to a two- Pomeron 
exchange of eikonal type. On the other hand, (|62|) was derived under the assumption that 
Yi and Y2 are large enough to justify the insertion of BFKL Pomerons between the loop and 
the impact factors, and we have to restrict AY to a region smaller than the total rapidity. 
If the total rapidity Y is much larger than the rapidity interval occupied by the loop - for 
example we could consider a loop of the finite size AY with AY ~ l/cjo ~ 4 - then such 
a "small" loops can be repeated many times and would play the role of the Pomeron self 
energy, leading to a renormalization of the Pomeron intercept. As mentioned before, the self 
energy is negative relative to the BFKL amplitude, and the renormalization therefore lowers 
the Pomeron intercept. If the absolute value of the renormalization is close or even larger 
than the 'bare' intercept uq = a{0) — 1, one may get close to the "critical" Pomeron or even 
obtain the "subcritical" Pomeron, as it has been discussed in ^, (see also |^ where 



a prescription for the renormalization of the supercritical (with a{0)p > 1) Pomeron was 
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proposed) . 

Before we do our estimate, it is useful to recall the nonperturbative renormalization 
caused by the pion loop insertion. For a single pion loop inside the Pomeron (at t = 0) we 
have a formula quite analogous to (62). Instead of Yi,Y2, convenient variables are: the size 
of the loop AY, the position of the center of the loop Yc = (Y — Yi + Y2)/2. Inside the pion 
loop we have the transverse momentum kt of the pion. Let us fix, for a moment, the value 
of Yc and calculate the pion loop contribution A'^: 

Here the first factor 3 results from the contributions of vr^, vr" and tt^. For the numerical 
estimate we neglect the pion mass , put the Pomeron pion coupling square gl = (77^,1 (sq) ~ 
20 mb and choose the mean transverse momentum of the first particle inside the Pomeron- 
pion vertex ko ~ 0.6 GeV. One needs this particle also to fix the rapidity of the vertex Yi 
(or Y2). The integral over Ay (|65[) is convergent, and for AY << Y the integration over Yc 
gives a factor Y. As a result, the one loop correction to the amplitude is equal to Y ■ A'^. 
Inserting two pion loops we obtain ^{Y ■ A'^)^, and summing over an arbitrary of loops we 
get the sum exp{Y ■ A'^) which means that the Pomeron intercept increases by A'^. 

Returning to the perturbative Pomeron loop insertion (62), we again choose = 0.3 
and assume Y >> AY. Dividing the amplitude (|6^ by the elastic forward amplitude (20) 
and considering only the integration over the loop size, we find: 

J [aAYY (47r)i7 ^ ' [AY^ 

Here we have cut the integration over the loop size at 6Y = AY. As the numerical coefficient 
here is extremely small, the only possibility to obtain a relatively large renormalization is 
to chose a very large AY >> l/cuo- On the other hand the loop renormalization A'°°^ is a 
NNLO BFKL effect and first we have to account for the NLO BFKL corrections which lowers 
the intercept down to tuo ~ 1/4. Therefore up to a very large AY the loop renormalization is 
still negligible; for ujq = 0.25 the values of A'°°^ and become comparable only for AY > 45 
(changing the value of ujq to the larger value ujq = 0.70 we are consistent with |T^, |2^). Of 



course, from the academic point of view we have to account for this renormalization effect, 
when s — > cxo; but at any reachable rapidity interval the value of A'°°^ is much less than ujq. 

Since the numerical value of the triple Pomeron coupling is not far from the nonper- 
turbative one, it is not surprising to see that the estimate ( |66D is also consistent with the 
nonperturbative evaluation. In the latter case we expect 

2 

A^7 ~ ^ J d{AY)dk^e''''^^ ^ 2 ■ 10-^5Ye''°^^ (67) 

for g^p ~ 1 GeV^i and k^ ~ 1 GeV^. 

Finally we mention that the estimate in (66) may still be a bit too large. Namely, recall 
that in our calculation of the Pomeron self energy and in our numerical estimate of its 
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magnitude we have restricted ourselves to the forward direction, that we have considered 
the values 1^ = I'a = 0. As discussed before, this is the point where the perturbative nature 
of the BFKL approximation becomes most visible, i.e. the 'distance' between pQCD and 
nonperturbative QCD is the largest. The mathematical manifestation is the l/u singularity, 
which immediately disappears if we depart from Zq, = 0. Therefore, as for the discussion of 
the numerical value of the triple Pomeron coupling, we expect that also our estimate of the 
renormalization due to the selfenergy may be slightly overestimated. 

6 Conclusions 

In this paper we have performed an estimate of the perturbative triple Pomeron vertex. 
Starting from the results of a leading-ln s analysis of QCD perturbation theory in the triple 
Regge limit, we have used the large- Nc limit to derive a fairly simple expression for the 
triple Regge inclusive cross section, which can be compared with the standard formulae 
used in the analysis of experimental data. A numerical estimate of the perturbative triple 
Pomeron coupling - which has a considerable theoretical uncertainty - indicates that its value 
is of the same order of magnitude as the nonperturbative one, obtained from earlier fits to 
experimental data. 

We have also tried to estimate the renormalization of the BFKL intercept due the 
Pomeron self energy loop. Formally speaking, this is NNLO effect and lies beyond the 
NLO corrections to the BFKL kernel. A numerical estimate - again with theoretical un- 
certainties - indicates that these corrections due to self interactions of the BFKL Pomeron 
are much smaller than the nonperturbative contribution (|65D . This is caused by the fact 
that this correction is proportional to the perturbative triple Pomeron vertex square, and 
it agrees with the evaluation of the Pomeron loop insertion based on the phenomenological 
value ( |63l) of the triple Pomeron vertex.. 

So finally we conclude that in spite of a 'huge number' fl = 7767 (see eq.(|5^) the 
perturbative triple Pomeron coupling is not large; it is in approximate agreement with the 
old phenomenological evaluations. 

As a future step, it would be interesting to understand better why numerical studies of 
the nonlinear evolution equations seem to find rather rapid saturation effects in spite of a 
rather small value of the perturbative triple Pomeron coupling. 

Appendix A: Counting factors of 2 and n 

In this appendix we give a brief summary of the normalization of the impact factors and 
the triple Pomeron vertex in perturbative QCD. Our starting point is eq.(4) which defines 
the impact factor. To be definite we consider the elastic scattering of two quarks (averaged 
over color and helicity of the incoming quarks). The lowest order diagram with color singlet 
exchange has two gluons in the t-channel (box diagram and its crossed counterpart), and in 
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the high energy hmit on finds 



Comparison with (4) yields the quark impact factor 



*'=2^*'^4]vr- '''' 

Equivalently, we could have defined the impact factor through the requirement that the 
leading order energy discontinuity should have the form 

r <Pk 1 

Next we turn to the triple Regge cross section. Eq.(23) defines the M^-discontinuity 
of a six-point function, a^. Turning again to lowest order QCD diagrams, we look at the 

diffractive process q + q ^ {qg) + qa, in the triple Rcggc limit; a typical QCD diagram is 
shown in Fig. 5. We define the triple Pomeron vertex through the LO equation 

^ f (Pq (Pk-i (fks^ 1 I, , ... , N 

I4{i-k,y^'l4{i + ksy^'- ^^^^ 

The analysis of the high energy behavior of the relevant QCD diagrams leads to the result: 



27V, ^ y (27r)3(q2)2^^9' '^'^'^ - k.f kl{l + k.f 



where V stands for the BFKL-type 2 — > 4 gluon vertex: 

V = (g')' t (73) 

{q-k,ni + k,-qf {q-k,fkl {I + k,y{l + k, ~ qf ^ ' 

Inserting the result for the quark impact factor we obtain the following normalization of the 
triple Pomeron vertex V: 

^^iW^^y (74) 

Appendix B: Impact factors in the conformal approxi- 
mation 

We discuss here briefly, with simple arguments and some approximations, the behaviour of 
the impact factors in conformal representation in the forward and non forward direction, but 
without giving a full momentum and conformal weight dependence. 
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Let us consider already the situation of zero conformal spin and work in coordinate 
representation. One can obtain the same results on studying the limiting case of the BFKL 
Pomeron eigenstate in momentum representation. 

We start from an impact factor which have some dominant support in a bounded region 
of size of order R in the coordinate space. We are interested in studying the behaviour of 
the following expression 

We are mainly interested in studying the behaviour in the small \l\ region. Therefore the 
main contribution in (|75|) will come from the integration in the region of large |ro|. It is 
therefore convenient to split the integration region according to |ro| < -R and |ro| > R. We 
shall be interested therefore in momenta \l\ < 1/R and neglect the first contribution. We 
note also the in the region ro| > -R it is a good approximation to consider |rjo| ~ l^ol for 
i = 1,2 since the external integral has support roughly for |rj| < R/2. We can therefore 
write, in a factorized form, 




2iu\ 



. (27r)2(27r)2' ^ ^' J\r,\>R {2'Kf 

and study the I dependence in g{v, I). The Tq integration gives 

q{vA) = nR 

' ' 2iu 

One can see that in the forward direction 1 = the second term does not give contribution 
and therefore g{v-,l) ~ l/{2iv). This leads infact to the correct behaviour of the BFKL 
Pomeron Green's function in forward direction, where under the v integration the integrand 
is not proportional to z/^, as, instead, in the non forward case. To analyze the limit of small 
IR? — >■ 0) one can keep the term (i?^Z^)^*^/(2ii/) and estimate with the saddle point method 
its contribution in such a limit. Again one can easily check that such a contribution is 
suppressed. For Z 7^ there is instead a cancellation of the v pole in the origin. Therefore 
we can write 

$(z/,0) ^ —^OF 

<l>(z/,Z) ^$o7VF, s>Z2>0 (78) 

We shall be interested typically, for the non forward case, to values of -R|Z| ~ 1, at the 
border of the approximations taken above to show the behaviour in (^). 
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